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Abstract 

Cyclidic nets are introduced as discrete analogs of curvature line 
parametrized surfaces and orthogonal coordinate systems. A 2-dimen- 
sional cyclidic net is a piecewise smooth C 1 -surface built from surface 
patches of Dupin cyclides, each patch being bounded by curvature 
lines of the supporting cyclide. An explicit description of cyclidic nets 
is given and their relation to the established discretizations of curva- 
ture line parametrized surfaces as circular, conical and principal con- 
tact element nets is explained. We introduce 3-dimcnsional cyclidic 
nets as discrete analogs of triply-orthogonal coordinate systems and 
investigate them in detail. Our considerations are based on the Lie 
geometric description of Dupin cyclides. Explicit formulas are derived 
and implemented in a computer program. 

1 Introduction 

Discrete differential geometry aims at the development of discrete equiva- 
lents of notions and methods of classical differential geometry. The present 
paper deals with cyclidic nets, which appear as discrete analogs of curvature 
line parametrized surfaces and orthogonal coordinate systems. 

A 2-dimensional cyclidic net in IR 3 is a piecewise smooth C 1 -surface as 
shown in Fig. Such a net is built from cyclidic patches. The latter 

are surface patches of Dupin cyclides cut out along curvature lines. Dupin 
cyclides themselves are surfaces in M 3 , which are characterized by the prop- 
erty that their curvature lines are circles. So curvature lines of the cyclidic 
patches in a 2D cyclidic net constitute a net of C 1 -curves composed of cir- 
cular arcs, which can be seen as curvature lines of the cyclidic net. 

'Partially supported by the DFG Research Unit "Polyhedral Surfaces" and the DFG 
Research Center Matheon 
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Discretization of curvature line parametrized surfaces and orthogonal co- 
ordinate systems is an important topic of current research in discrete differ- 
ential geometry. The most recognized discretizations (see |BS08j) include: 

• circular nets, that are quadrilateral nets with circular quadrilaterals, 

• conical nets, that are quadrilateral nets, such that four adjacent quadri- 
laterals touch a common cone of revolution (or a sphere) , 

• principal contact element nets, that are nets built from planes contain- 
ing a distinguished vertex (such a pair being called a contact element), 
such that the vertices constitute a circular net and the planes form a 
conical net. 

The most elaborated discretization of curvature line parametrized surfaces 
are circular nets. They appear indirectly and probably for the first time in 
|Mar83j . The circular discretization of triply orthogonal coordinate systems 
was suggested in 1996 and published in Bol>99 . The next crucial step in the 
development of the theory was made in |CDS97| . where circular nets were 
considered as a reduction of quadrilateral nets and were generalized to arbi- 
trary dimension. An analytic description of circular nets in terms of discrete 
Darboux systems was given in |KS98j . and a Clifford algebra description of 
circular nets can be found in [BHJOlJ. The latter was then used to prove 
the C°°-convergence of circular nets to the corresponding smooth curvature 
line parametrized surfaces (and orthogonal coordinate systems) in [BMS03J . 
Conical nets were introduced in |LPW + 0o] as quadrilateral nets admitting 
parallel face offsets. This property was interpreted as a discrete analog of 
curvature line parametrization. Circular and conical nets were unified un- 
der the notion of principal contact element nets in }BS07] (see also [P W07] ) . 
The latter belong to Lie sphere geometry. 

All these previous discretizations are contained in cyclidic nets. The ver- 
tices of a cyclidic patch are concircular, therefore the vertices of a cyclidic 
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net build a circular net. On the other hand the tangent planes at these ver- 
tices form a conical net, so vertices and tangent planes together constitute 
a principal contact element net. A cyclidic net contains more information 
than its corresponding principal contact element net: The latter can be seen 
as a circular net with normals (to the tangent planes) at vertices. At each 
vertex of a cyclidic net one has not only this normal, but a whole orthonor- 
mal frame. The two additional vectors are tangent to the boundaries of the 
adjacent patches (cf. Fig. 1.1), resp. the patches are determined by this 
condition. We explain in detail how a cyclidic patch is determined by its 
vertices and one frame, and derive its curvature line parametrization from 
that data. A cyclidic net is uniquely determined by the corresponding circu- 
lar net and the frame at one vertex. It is worth to mention that this frame 
description played an essential role in the convergence proof in [BMS03J. 

Composite C 1 -surfaces built from cyclidic patches have been considered 
in Computer Aided Geometric Design (CAGD) |Mar83l IMcL851 IMdPS861 
IDMP931 ISKD96| . A large part of the book [JNM88] by Nutbourne and 
Martin deals with the Euclidean description of "principal patches" and in 
particular cyclidic patches. It also contains parametrization formulas for 
cyclidic patches, determined by a frame at one of the four concircular ver- 
tices. Unfortunately these formulas are rather complicated and turned out 
to be not appropriate for our purposes. A sequel, focused mainly on surface 
composition using cyclidic patches, was planned. However, the second book 
never appeared, probably due to the limited possibilities of shape control. 
Another application of Dupin cyclides in CAGD is blending between dif- 
ferent surfaces like cones and cylinders, which are special cases of Dupin 
cyclides (see for example [Deg02 j , which also contains a short analysis of 
Dupin cyclides in terms of Lie geometry) . 

In contrast to the previous works on surfaces composed of cyclidic patches, 
the present work relies essentially on Lie geometry. We use the elegant 
classical description of Dupin cyclides by Klein and Blaschkc [Klc26, Bla29j 
to prove statements about cyclidic patches, and to derive their curvature 
line parametrization. On the other hand, cyclidic nets as well as circular 
nets are objects of Mobius geometry, and we argue several times from a 
Mobius geometric viewpoint. In Appendix [A] we give a brief overview of the 
projective model of Lie geometry. 

Another central part of the present work is a new discretization of or- 
thogonal coordinate systems as cyclidic nets (in particular in M 3 ) . Smooth 
orthogonal coordinate systems are important examples of integrable systems 
[Zak98] . In the discrete setting analytic methods of the theory of integrable 
systems have been applied already to circular nets in [DMS98J (5-method) 
and in |AVK99| (algebro-geometric solutions). Our approach to the dis- 
cretization is the following: In the smooth case the classical Dupin theorem 
states that in a triply orthogonal coordinate system of M 3 all coordinate 
surfaces intersect along common curvature lines. This translates to a 3D 
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cyclidic net with its 2D layers seen as discrete coordinate surfaces, for which 
orthogonal intersection at common vertices is required. It turns out that 
this implies orthogonal intersection of the cyclidic coordinate surfaces along 
common discrete coordinate lines (cf. Fig. 1.2). Moreover it is shown that 
generically the discrete families of coordinate surfaces can be extended to 
piecewise smooth families of 2D cyclidic nets, which gives orthogonal coor- 
dinates for an open subset of M 3 . The corresponding theory is extended to 
higher dimensions. 




Figure 1.2: 2D layers of a 3D cyclidic net intersect orthogonally in discrete 
coordinate lines, and in particular at vertices. 

Our construction of cyclidic nets, based on the projective model of Lie 
geometry, has been implemented in a Java webstartQ All 3D pictures of 
cyclidic nets have been created with this program. A brief description of 



the implementation is given in Section 3.5 



Notations. We denote the compactified Euclidean space 

R N := R N U {oo} S N . 

Homogeneous coordinates of a projective space are also marked with a hat, 
i.e. x = [x] = P(Kc). The span of projective subspaces Ui = P(J7j),i = 
1 , . . . , n we write as 

inept, ...,U n }:= P(span(f7i, . . . , U n )). 

The corresponding polar subspace with respect to a quadric is denoted 

pol[l7i, ...,[/„]:= pol[inc[C7i, ...,U n ]} = P(span(ft, . . . , f/„) ± ), 

where _L denotes the orthogonal complement with respect to the inner prod- 
uct defining the quadric. 



Available at www. math, tu-berlin/geometrie/ps/sof tware . shtml 
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Pay attention to the different use of (■,■) in our work. We denote this 
way the inner product on M N+1,2 of signature (N + 1,2) which defines the 
Lie quadric in the projective space MP^" 1 " 1,2 = P(M 7V+1 ' 2 ) (cf. Appendix [a| , 
as well as the standard Euclidean scalar product in K . The products can 
be distinguished by the notation of the involved vectors: As homogeneous 
coordinates of MP JV+1 ' 2 the vectors in ~R N+1 ' 2 carry a hat, while vectors in 
M. N don't. Thus (xi,X2) denotes the product in ]R Ar+1 ' 2 and (x\,X2} denotes 
the Euclidean scalar product in R . 



2 Dupin cyclides 

Dupin cyclides are surfaces which are characterized by the property that all 
curvature lines are generalized circles, i.e. proper circles or straight lines. 
They are invariant under Lie transformations, and in this paper we use the 
projective model of Lie geometry for their analytic description [Kle26, Bla29, 
IPin85j . 



2.1 Description in the projective model of Lie geometry 

In the classical book [Bla29] by Blaschke one finds 

Definition 2.1 (Dupin cyclide). A Dupin cyclide in M 3 is the simultaneous 
envelope of two 1 -parameter families of spheres, where each sphere out of 
one family is touched by each sphere out of the complementary family. 

From the Euclidean point of view there are 6 different types of Dupin cy- 
clides. A Dupin cyclide is called to be of ring, horn, or spindle type depend- 
ing on whether it has 0, 1 or 2 singular points. Moreover one distinguishes 
bounded and unbounded cyclides, where the latter are called parabolic. The 



standard tori are Dupin cyclides (cf. Fig. 2.1 ) and each other cyclide can be 
obtained from a standard torus by inversion in a sphere [Pin86] . In particu- 
lar the parabolic cyclides are obtained if the center of the sphere is a point 
of the cyclide (a parabolic horn cyclide with its singular point at infinity is a 
circular cylinder, a parabolic spindle cyclide with one of its singular points 
at infinity is a circular cone). 






Figure 2.1: The standard tori are special types of Dupin cyclides. 

The class of Dupin cyclides is preserved by Mobius transformations. In 
the Mobius geometric setting only the three types "ring", "spindle", and 
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"horn" remain, as there is no distinguished point at infinity. So far, this can 
be found in the compact survey [Pin86j by Pinkall. 

In Lie geometry, which contains Mobius geometry as a subgeometry, the 
distinction between the different types of cyclides, where the type is deter- 
mined by the number of singular points, becomes obsolete: There are Lie 
sphere transformations which map points to proper spheres and vice versa, 
in particular the number of singular points of a Dupin cyclide is not a Lie 
invariant. The following elegant characterization of Dupin cyclides is a re- 



formulation of the Definition 2.1 in the projective model of Lie geometry. 
It can be found implicitly already in [Kle26] and in more detail in |Bla29j . 
while in |Pin85] the analogous statement is proven for the general case of 
so-called Dupin hyper surf aces. 

Theorem 2.2. A Dupin cyclide in R 3 corresponds to a polar decomposition 
o/IRP 4 ' 2 into projective planes P^ and P^ of signature (+ H — )• 

Indeed oriented spheres in Lie geometry are described as points on the 
Lie quadric 

Q 4 ' 2 = P(L 4 ' 2 ) C MP 4 ' 2 where L 4,2 = {v £ M 4 ' 2 | (v, v) = 0} 

and spheres si,S2 in oriented contact are modelled by polar points, i.e. 
Si = [vi] with (vi,V2) = (cf. App endix |A|). Thus talking about polar 
planes P^ and P^ in Theorem 2.2 is another description of two families 



of spheres as in Definition |2.1| This motivates 
Definition 2.3 (Cyclidic families of spheres). Denote 

V^.y C = \Planes in MP 4,2 with signature (+ H — )} • 
A 1 -parameter family of spheres corresponding to a conic s ectio n P n Q 



■1.2 



with P S Vtyc is called a cyclidic family of spheres (cf. Fig. 2.2). 



Two spheres s\, S2 in oriented contact span a contact element, which con- 
sists of all spheres in oriented contact with s± and S2 and is modelled as 
the projective line inc[si,S2] contained in the Lie quadric Q 4 ' 2 . Such lines 
L C Q 4 ' 2 are called isotropic lines and we define 

J*?q' 2 = {isotropic lines in MP 4,2 } = |contact elements of M 3 | . 

To be precise, the Lie- geometric char acte rization of Dupin cyclides as ele- 
ments of Vtyc coincides with Definition 2.1 if one considers oriented spheres 
in oriented contact and reads "spheres" as "generalized oriented spheres" 
(i.e. hyperspheres, hyperplanes, and points). A circle then is also a Dupin 
cyclide: One family of spheres consists of points (i.e. spheres of radius 0) 
constituting the circle, the other family consists of all oriented spheres con- 
taining this circle (including the two oriented planes). This coincides with 
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Figure 2.2: Polar planes p( 1 ) i p( 2 ) g vtyc describe a Dupin cyclide C: Each 
plane intersects the Lie quadric Q 4,2 in a non-degenerate conic C^' = Q 4 ' 2 D 
?W , and each conic corresponds to one of the 1-parameter families of spheres 
enveloping C. Each point x of the cyclide is obtained as contact point of 
spheres E C^ 1 ' and s^ 2 ) E C^ 2 \ the contact elements of C are given as 
inc[sM, sOO] c Q 4 ' 2 . 



Lemma 



A.4 



which states that points in M 3 are concircular if and only if 
their representatives in the Lie quadric are coplanar. Note the implication: 
if three spheres in a cyclidic family are point spheres, i.e. spheres of van- 
ishing radius, then the whole cyclide has to be a circle. This corresponds 
to the fact that a Dupin cyclide in Mobius geometry, which is a 2-surface, 
contains at most two singular points. If one treats circles as Dupin cyclides, 
the correspondence described in Theorem 2.2 is indeed a bijection. 

To simplify the presentation we exclude the circle case in the following 
considerations, which allows us to see Dupin cyclides as surfaces in 3-space 
with at most two singular points. 



2.2 Curvature line parametrization of Dupin cyclides 

Proposition 2.4. Curvature lines of Dupin cyclides are circles, each such 
circle being contained in a unique enveloping sphere. In particular the en- 
veloping spheres of a Dupin cyclide are its principal curvature spheres. 

Given a cyclide C and its enveloping spheres CW = pW n Q 4,2 , i = 1,2, 
points along the curvature line on s" E C^ 1 ' are obtained as contact points 



of sW with all spheres of the complementary family C^') (cf. Fig. 2.2). 
Curvature lines from different families intersect in exactly one point. 

Proof. Let gW E be a fixed sphere which is not a point sphere. By 



Theorem 2.2 all contact elements of C containing s^' are of the form 
inc[sW, s^], s& E i ^ j. They describe a circle on sW since 
is contained in the plane P^ (in this case inc is 



3-dimensional and one can use Proposition A. 5). Such circles are curvature 
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lines, since if a surface touches a sphere in a curve, this curve is a curvature 
line and the fixed sphere is the corresponding principal curvature sphere. 
Two curvature circles C s''' G CW, i = 1,2, intersect in the unique con- 
tact point of and s^ 2 \ Since each non-singular point of the cyclide is 
contained in a unique contact element, all curvature lines are circles. □ 



Proposition 2.4 implies that parametrizations of polar cyclidic families 
of spheres induce a curvature line parametrization of the enveloped Dupin 
cyclide (cf. Fig. fllty . 



Lemma 2.5 (Parametrization of conies). Let si,S2,S3 be three points on a 
non- degenerate conic C in a projective plane P and let 

I a 12 ai 3 N 
<3 = ai2 a 23 
\ai 3 a 2 3 

be the matrix representation of the corresponding quadratic form q in the 
basis (si, ^2> S3) ; ^.e. s i = [1>0,0] T , s 2 = [0, l,0] T ,s 3 = [0,0, 1] T and aij = 
q(si,'sj). Then s = [s] : R — >■ P defined by 

( a 23 • (2t 2 - 3t + 1) \ 
s(t)= a 13 -(t-t 2 ) (2.1) 
V ai2 • (2t 2 -t) ) 

is a parametrization of C for which s(0) = si, s~(i) = s 2 , s(l) = S3. 

This is a standard result obtained by projecting a conic C from a point 
p £ C to a line, see [BEG02] for example. A direct computation shows 
s T Qs = for all t 6 M, as well as the proposed normalization. 

Applied to cyclidic families of spheres this yields the following Proposition 
(according to later application it is convenient to use upper indices i): 

Proposition 2.6 (Parametrization of a cyclidic family of spheres). Let 

q(i) — p(i) p| Q4,2^p(«) g Vcyc, be a cyclidic family of spheres which con- 
tains spheres s\ , s% 9, s 3 i . A parametrization = [s®] '■ K — > which 
maps 1 — y s\ , \ \- > s| , 1 1— > s| , is given by quadratic polynomials 



^(t) = (sf\si\si ) ). 



( 45 • (2t 2 - 3t + 1) 



*23 

g • " * 2 

V ag-(2t 2 -t) 



(t-t 2 ) • (2-2) 



.ffere the columns sfp are normalized homogeneous coordinates (A.l) and 



a jk ■— V s j ' s fc /• 
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Proof. The parametrization (2.2) follows immediately from (2.1) since the 
points s% \ Sg^ 6 span PW. (There are no three collinear points 
on CW due to the signature (+ H — ) of PW, which means that CW is a 
non-degenerate conic). □ 

Theorem 2.7 (Curvature line parametrization of Dupin cyclides). Let 
S W — [jW] ; k — > C^\i = 1,2, oe parametrizations (2.2) of the cyclidic 
families of spheres enveloping a Dupin cyclide C C K 3 . A (Lie geometric) 
curvature line parametrization of C in terms of contact elements is given by 



L : 



1,2 



(tLtaJ^inc^Cti).^^)]. 



(2.3) 



The corresponding Euclidean curvature line parametrization is given by the 
first 3 components of the induced parametrization 



L 



1,2 



(tl,t 2 ) ^ 



rWfa^it^-rWjh^ih) 
r( 2 )(t 2 )-rM(ii) 



(2.4) 



The right hand side of (2.4) are normalized homogeneous coordinates (A. 2) 



of the unique contact point of ( |2.3[ ) , which is the corresponding point of the 
cyclide. The oriented radius r^(U) of s®(ti) is given by 



„(i) 



+ 
+ 



l 23 
l 31 
l 12 



(2if - 3ti + 1) 
(2tf-i), 



(2.5) 



where (s^) r 



is £/te e r -component o/s^- and = (s'j\'s'^) 



Remark. Note that x can be constant along parameter lines (corresponding 
to singular points of a cyclide), but L cannot. 



Proof. Definition 2.1 translates to (2.3), and because of Proposition 



is a curvature line parametrization of C in Lie geometry. Lemma A. 3 



2.4 this 



(2.4), and a direct computation using (2.2) leads to (2.5) 



implies 
□ 



2.3 Cyclidic patches 

Definition 2.8 (Cyclidic patch). A cyclidic patch is an oriented surface 



patch, obtained by restricting a curvature line parametrization (2.4) of a 



Dupin cyclide to a closed rectangle L\ x I2. We call a patch non-singular if 
it does not contain singular points. 

Geometrically, a cyclidic patch is a piece cut out of a Dupin cyclide along 
curvature lines, i.e. along circular arcs, as in Fig. |2.3| If the vertices of a 
patch are non-singular points, the boundary consists of four circular arcs 
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which intersect orthogonally at vertices. Moreover, a patch is non-singular 
if and only if its boundary curves intersect only at vertices. 

Now we introduce a notation which will be convenient for the description 
of cyclidic nets in Section [3J If / : [0, 1] x [0, 1] — > M 3 is a (parametrized) 
cyclidic patch, we denote its vertices by 

x = /(0, 0), xi = /(l, 0), xi 2 = /(l, 1), x 2 = /(0, 1). 

The notation (x, x\, x\ 2 , x 2 ) refers to the vertices in this order. Further we 
write 5 iX — X 1 X and S iX j — X^j X j respectively. Finally we denote by 
x^x~i the boundary arc of a cyclidic patch connecting the vertices x and X{. 

Definition 2.9 (Vertex frames and boundary spheres of a cyclidic patch). 
Let f be a cyclidic patch with vertices (x, x±, X12, x 2 ). The vertex frame of f 
at x is the orthonormal 3-frame B = (t^\ t^ 2 \ n), where £W is the tangent 
vector of x^Xi directed from x to xi and n is the normal of the supporting 
cyclide. The boundary spheres of f are the principal curvature spheres 
supporting the boundary arcs. 

The vertex frames always adapt to the boundary arcs of a patch. For 
example the frame B\ = (t^\t^\ni) at x\ contains the unit vectors t± 
and t^\ where tj is the tangent vector of x^x\ directed from x\ to x, and 



t\ is the tangent vector of x\,xi 2 directed from x\ to x\ 2 (cf. Fig. 2.5) 



The perpendicular bisecting hyperplane of the segment [x, Xi] we denote 
h {i) = [ P eR N I (p,5 iX ) = (x + h iX ,5 t x) = ^(|M| 2 - ||x|| 2 )}. (2.6) 

By IJW we denote the reflection in fcw which interchanges x and x%. A 
vector v at x is accordingly mapped to V{ at due to 

(7) (7) 

The corresponding reflections in the hyperplanes h- are denoted by Hf . 



(cf. Fig. 2.4). These mappings extend naturally to tuples of vectors and we 



denote the image of a frame B by H^(B). 
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Lemma 2.10. Let x, xi,x±2, X2 be four distinct points on a circle. Then the 



reflections (2.7) in hyperplanes (2.6) satisfy 



(2.8) 



Proof. Since the considered points are concircular they lie in a 2-plane II. 
All hyperplanes intersect in the affine orthogonal complement of II through 
the center of the circle (cf. Fig. 2.4). □ 




Figure 2.4: Perpendicular bisecting hyperplanes between concircular points 
intersect in the axis of the circle. 



Proposition 2.11 (Geometry of cyclidic patches). For a cyclidic patch with 
vertices (z, x\, xi%, x%) one has: 

i) The vertices of the patch lie on a circle. 

ii) Vertex frames at neighboring vertices are related by reflection in the bi- 
secting plane (2.6), possibly composed with the direction reversion of the 
non- corresponding tangent vector. The change of orientation depends 
on whether the vertex quadrilateral (x, xi, X12, Xz) is embedded or not. 
For example frames B and Bi at vertices x and X{ are related as follows: 

• If the points x and Xi on the circle are not separated by Xj and 
x\2, then 

Bi = H®(B). (2.9) 

• If x and Xi are separated by xj and x\2, then 

Bi = (H^ o F®)(B) = (F (j) o H®)(B), i + j (2.10) 
where F^' denotes the orientation change of in {t^ l \t^ 2 \n). 

In particular lines x + Mn and x% + Mn, intersect in the center of the 
boundary sphere containing the boundary arc x^Xi (cf. Fig. 2.5). 
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Figure 2.5: Geometry of a cyclidic patch. 



Proof, i) Consider the cyclidic families of spheres C^\C^ 2 ' enveloping the 
cyclide supporting a cyclidic patch /. Denote the boundary spheres of / 



by s^jS^ £ C^ 1 ) and s^ 2 ^,s^ 6 C^ 2 ). The contact elements at vertices of 
/ are given by L = incjs^ 2 ), sW], L\ = inc[ 



L 2 = vac[a$\aV>] (cf. Fig. |2.6|. The space V 



mc|sW,4 2) ],i>l2 



(i) 



, (2) (1) 

mc[L,Li,Li2,L 2 ] C 

Since V contains 



and 

4,2 



is at most 3- dimensional, as V = mc[s 
proper spheres, the intersection V D P(e^~) is at most 2-dimensional. 
Lemma A. 4 the vertices (x, x\, x%2, £2) 



By 



are concircular. 



>12 




Figure 2.6: Combinatorics and geometry of contact elements (as isotropic 
lines in Q 4,2 ) at vertices of a cyclidic patch. 

ii) The sphere is symmetric with respect to which implies (n) = 
m. Moreover H^(t^) = tf\ since each circular arc containing x and X{ 
is symmetric with respect to h®. As is an isometry, orthogonality is 
preserved and H^(t^) = rti^ for i 7^ j. Thus i? and i?« are necessarily 



related by (2.9) or (2.10) and it remains to distinguish those cases: Since we 



know that a cyclidic patch contains either 0,1, or 2 singular points, Figs. 2.7 



and 2.8 show all possible cases (up to Mobius transformations). It turns out 



that H®(tW) = i.e. pT0| ) holds, if and only if B and B, t are of same 

orientation. This happens if and only if 0C . X 2 ell e separated by Xj, x\i on the 
circle through vertices. □ 



Spherical patches. A sphere s may be seen as the limit e — > of a 
family C e of Dupin cyclides. In this case exactly one of the t wo c orresponding 
families of conies C Q 4 ' 2 converges to s G Q 4 ' 2 (cf. Fig. 2.2). We call the 
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Figure 2.7: If the vertices (x, x\,x±2, £2) of a cyclidic patch build an embed- 
ded quadrilateral, vertex frames at adjacent vertices have opposite orienta- 
tion. 




Figure 2.8: Vertex frames of a cyclidic patch, for which the vertices 
(x, x\, x±2, X2) build a non-embedded quadrilateral. 



corresponding limits of cyclidic patches spherical patches and consider them 
as degenerated cyclidic patches for a unified treatment. Continuity implies 
that a spherical patch also has circular boundary arcs, and that the crucial 



geometric properties of cyclidic patches described in Proposition 2.11 hold as 



well (cf. Figs. 2.9 and 2.10). Note that the supporting sphere is determined 



by the boundary curves of a spherical patch, while in contrast the boundary 
curves of a generic cyclidic patch are never contained in a single sphere. 






Figure 2.9: Proposition 2.11 also holds for spherical patches. 



Proposition 2.12. For a spherical patch there is a unique pair of orthogonal 
1 -parameter families of circles on the supporting sphere, determined by the 
condition that the circular boundary arcs of the patch are contained in circles 
of those families. 



Proof. Consider a spherical patch / contained in a sphere, and the four 
circles determined by the boundary arcs of /. First suppose that the four 
circles are distinct and consider one pair of opposite circles. Those two circles 
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Figure 2.10: Combinatorics and geometry of contact elements (as isotropic 
lines in Q 4 ' 2 ) at vertices of a spherical patch. The Vertices are concircular, 
i.e. their representatives in the Lie quadric are coplanar, and the contact 
elements span a 3-space. 



share either 0, 1, or 2 points. These cases can be normalized by a Mobius 
transformation to the cases of two concentric circles, two parallel lines, or 
two intersecting lines. In each case the two families of orthogonal circles are 



uniquely determined (cf. Fig. 2.11). If two boundary circles coincide, the 
two others have to be distinct and one argues with the same normalization. 

□ 




Figure 2.11: There are two types of orthogonal pairs of 1-parameter families 
of circles. The type of a pair is determined by two circles from one family. 



A point x of a spherical patch is a singular point if all circular arcs of one 



family pass through x (cf. Fig. 2.11 ). Considering the two families of circles 



on the whole sphere, one has exactly 2 or 1 points of this type. 



Proposition 2.12| has the following consequence: 



Corollary 2.13. For a spherical patch there exist orthogonal parametriza- 
tions with all coordinate lines being circular arcs. The parameter lines are 
uniquely determined by the patch vertices and the boundary arcs. 
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Cyclidic patches determined by frames at vertices 



Theorem 2.14. Given four concircular points, there is a 3-parameter family 
of cyclidic patches with these vertices. Each such patch can be identified with 
an orthonormal 3-frame at one vertex. 



Proof. Due to Proposition |2.11 it is enough to show that for concircular 
points X = (x, xi, X12, X2) and any orthonormal 3-frame B = (t^\t^ 2 \n) 
there is a unique cyclidic patch with X as vertices and vertex frame B at x. 
Denote by cx the circle through X. First obtain frames Bi, B12, B<i at 



xi,xi2,X2 from B by successive application of (2.9) and possibly (2.10), 
depending on the ordering of the points X on cx- The frames are well 
defined: In the case that X are vertices of an embedded quad, this follows 
immediately from Lemma 2.10 In the non-embedded case the additional 



orientation changes involved in (2.10) cancel out (B, Bi are related by (2.10) 
if and only if Bj , B^ are as well) . 

The frames at vertices induce contact elements C = (L, Li, L12, L2), where 
(x,n), etc. D ue to (2.9), (2.10), contact elements L,L*i share a sphere 
The intersection point 



L = (x,n), etc. D ue to (|2.9|) 
as depicted in Fig. 



2.12 



i'*' = LOLj with unoriented radius r^) 



is the center of a sphere 
= |U,_ C W||. Note that 









f \y\ 




PS 




/ \ 










N 





Figure 2.12: Construction of proper boundary spheres of a cyclidic patch. 



(0 



x,Xi and B already determine a circular arc x^Xi C M 3 . Since x,xi G s 
and is tangential to at x, we have x~xi C sW. 

The generic case. Generically, i.e. if and only if x + Mn and the axis of 
cx are skew, C is as in Fig. 2.6. Fix i,j G {1, 2} such that i 7^ j and denote 
£;« the unique circle containing the circular arc £^~Xi determined by x, Xi 
and iW £ B. First we show that there is a unique cyclide containing the 
contact elements C, for which fe® C i s a curvature line: 
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Denote by V C pol[sW] t he 3- dimensional subspace of . 



p4,2 C orrespond- 
every cyclide 



2.2 



ing to kfr' (cf. Proposition A. 5). According to Theorem 
can be described by a plane G Pcj/c (determining pM via duality). The 
corresponding cyclid e to uches in feW if and only if V = inc[P^,sW], 
see also Proposition 2.4, and this determines P^ uniquely because of the 



following: Due to the signature (+H — ) 
not be in oriented contact. Therefore 



of pW, the spheres sW and can- 



G V implies V ^ pol[s^ ]. Since 
dimF = 3, dimPW = 2, and pW) c pol^ ], it follows = VTlpol^]. 
The obtained cyclide C is independent of the choice of i G {1,2}, since k^' is 
also a curvature line of C: We know that D fe«) is a curvature sphere of 
C and therefore C touches s^' in a unique circle k. This circle is tangent to 
in x due to the orthogonality of curvature lines. Thus Xj G C D s^' im- 
plies Xj G fc, and we have A; = A;^' as A;^' is the unique circle through x and 
which is tangent to i«) in x. According to this, all four circles containing 
the circular arcs determined by vertices X and frames (B, B\, B12, B2) are 
curvature line of C. 

Finally the frame B at x determines the orientation of curvature lines of 



C, and hence a unique cyclidic patch for the given data (cf. Fig. 2.13). 




Figure 2.13: Vertices and vertex frames of a generic cyclidic patch determine 
the supporting cyclide and the patch itself. 



The spherical case. The contact elements C are as in Fig. 2.10 and share 
exactly one common sphere if and only if x + In intersects the axis of the 
circle cx or is parallel to it. Clearly there is a 1-parameter family of such 
normals n, and they are in bijection with all oriented spheres containing cx- 
Denote the corresponding 2-parameter family of frames by Bs- As there is 
a 3-parameter family B of orthonormal frames at x, for a frame B G Bs it is 
possible to choose a continuous variation of generic frames B £ G B\Bs with 
B £ ^> B. According to our consideration of spherical pages (cf. pp. 12 1) 
this induces a continuous variation of cyclidic patches with fixed vertices 
and gives a unique spherical patch in the limit. □ 



Next we describe the restriction of cyclidic patches to subpatches, deter- 
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mined by the choice of points on the boundary arcs of the initial patch. 



Proposition 2.15. For a cyclidic patch with notation as in Fig. 2.1J\ 
points xi and x% on adjacent boundary arcs determine a unique cyclidic 
subpatch. If xi and X2 are non-singular points, the subpatch has four ver- 



tices (x,xi,xi2,X2) which are concircular, and its vertex frames satisfy (2.9) 



resp. (2.10) depending on the ordering of vertices. 




Figure 2.14: Choosing points on the boundary determines a subpatch. 

Proof. The generic case. Let / : [0, 1] x [0, 1] — > M. 3 be a curvature 
line parametrization of the patch for which x = /(0,0),Si = f(uo, 0),X2 = 
f(0,v ),0 < u ,v < 1. This determines a unique subpatch / = /|[ 0)Uo ] X [ 0) „ ]- 
The fact that x,x\,X2 are non-singular implies that x\2 '■= f(uo,vo) is also 
non-singular and in particular distinct from the initial three points. Hence 



the claim follows from Proposition 2.11 

The spherical case. This is proven analogous to the spherical case of 



Theorem 2.14, using a variation of frames B £ at x. Note that this time one 



has to consider an additional variation of points (xi) £ inducing a variation 
f £ of generic subpatches, for which only the common vertex x is fixed. □ 



Curvature line parametrization of cyclidic patches 



According to Theorem 2.14, circular vertices X = [x,x\,x\2,X2) and an or- 
thonormal 3-frame B = (t^ , t^ , n) at x determine a unique cyclidic patch. 
In the generic case one obtains an Euclidean curvature line parametrization 
as follows: 



0) Identify the notation used in the context of Dupin cyclides, with the 
notation for cyclidic patches: Boundary spheres s^,sf ,i ^ j of a 



(0 „(*) 



used in Proposition 



2.6 



The 



patch correspond to spheres s\ 

(i) 

third sphere s^ needed for a parametrization (2.2 ) of a cyclidic family 
of spheres is now denoted . 



(2) a (2) ,( 2 h 



1) Determine Lie coordinates (A.l) of the spheres (s' ',<t' >,s 



and 



scribed in Proposition 2.16 



from the data (X, B). How this can be done is de- 
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2) Use the coordinates obtained in step 1) in (2.2 ). This gives parametriza- 
tions s^ 1 - 1 , s^ 2 ) of the cyclidic families of spheres enveloping the support- 
ing cyclide. 

3) Putting parametrizations s^ 1 ^ and s^ 2 ) from step 2) in (2.4) yields a 
parametrization x of the whole cyclide. 

4) A Euclidean curvature line parametrization of the patch is given by 
the first three components of the restriction of x to [0, 1] x [0, 1]. In 
particular (|, |) is mapped to the contact point of and a^ 2 \ 



All 3D pictures in this work except Fig. 2.1 have been produced with 
a program which implements the above described algorithm (see also Sec- 
tion 



3.5). 



Proposition 2.16 (Lie coordinates for curvature spheres of a cyclidic patch). 
Let f be a generic cyclidic patch with concircular vertices (x, x%, x%2, £2) cmd 
vertex frame B = (tP~\t^ 2 \n) at x. Further denote 5 %X — X% X • 

If n JL Six, Lie coordinates (A.l) of the boundary sphere «W f f are 



3(0 



x + 
+ 



\Six\ 



2(5iX,n] 
Ifccll 2 



-n + e + 



\x\\ 2 (5iX,n) + ||5jx|| 2 (x, n) 
(5iX,n) 



+ 



2(5iX, n) 



(2.11) 



Lie coordinates (2.11) of the shifted boundary sphere 8j, i 7^ j, are obtained 
using the shifted formula: (5ix) = x\2 — Xj, n — V nj, where 

the shifted normal Uj at Xj is given by 



2 {6jX,n) 5.x. 
(5jX,5jx) 



(2.12) 



If t® if^i^, Lie coordinates (A.2) of the midpoint y® of the boundary arc 
x^Xi C s" are given by 



where 



,(') 



x + p®v® + e + \\x + p yv v 



(0„W||2 f 



|fe||t (i) + 5iX 
UixW 2 



I fell 



(2.13) 

(2.14) 
(2.15) 



r 2(5iX,vW) 2(6ix,t®) +2\\Six\\' 

For the enveloping sphere erCO of f which touches s(0 in yW one obtains 



a 



{^ffl^-(s^ffly^ 



(2.16) 



Generically a^i) is a proper sphere and its Lie coordinates ( A.l ) are obtained 
by normalizing the eg- component of (2.16) to 1. 
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Proof. Denote h^' the bisecting plane (2.6). If n JL 5{X, the cent er 
_l_ r M n Q f S M is the intersection x + RnflM'' as shown in Fig. 



2.12 



Solving 



-. Since llrel 



1. 



(x + r^ra, 6{x) = (x + \5iX, 5{x) for rW gives n 1 ' _ 2 ^ — ^ 

the value is indeed the radius of s" and one easily checks that positive 
radius corresponds to the inward normal field. Putting c^' and in (A.l) 



gives (2.11). Here we used 



,(0||2 



x|| 2 + 2r (i) (x,n) 



|j;|| 2 (5j3;, n) + ||<5j:r|| 2 (2;, n) 
(5iX,n) 



The shifted normal at Xj is obtained by reflection of n in the bisecting 
plane between x and Xj, i.e. Km holds (cf. (pTn)). 

If £W || the midpoint yM of 3^ ^ X<i IS the i nterse ction of the angular 
bisector between 5iX and with M'' (cf. Fig. 2.15-left): Let a G (0, 7r) 













/ a r 


c 


?C\ r 





7 = 




Figure 2.15: Construction of the midpoint of a boundary arc. 



be the angle between £W and <5jX, and denote c the center of the circle 
determined by x,Xi and its tangent x + MtW. Then the angle Z(:r,c, yW) 

■right for a G (0,^) and Fig. 2.15 -middle for 
c||, which implies j3 



also equals a (cf. Fig. 
a G (|,7r)). Clearly ||x 



2.15 



and thus 



7 = f - p = § . So y« = x + Ryd) n fcW with V W = \\6ix\\t® + 5 iX as 
pW|| = 1. This time solving {x + p^v^,5ix) = {x + ^5iX,5ix) for /jW leads 
to 

ll(5,-x|| 2 ll&cll 



(0 



(0 



2(5iX,vW) 2(6iX,t^) + 2\\5iX 
Thus for fw [f5jX the Lie coordinates (A. 2) of yW are exactly (2.13). 



The last formula (2.16) follows immediately from ay' 

(cf. Fig. 



mc s 



(0 



,y w ]n 



2.16). The solution is 



pol[s^] by solving (Xs® + /iyW,s^- 

unique (up to normalization of homogeneous coordinates) as and are 
not in oriented contact for generic patches (cf. p. 35). □ 
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mc s 




c pol[«W] n pol[s 



Wi 




,(0 



Figure 2.16: There is a unique sphere o - ") touching in 



which is also 

. This determines the plane 6 Pqyc- (Note 
that the shown configuration is 4-dimensional) . 



in oriented contact with s 



Remarks. 1. Note that (2.16) can be applied using any point from 



the interior of the boundary arc x, a^. 



#0 



If n _L the sphere becomes a plane with coordinates s 

, i^eoo+e?- . This is obtained for example by renormalizing (12.111) 



via multiplication with 



2(&iX,n) 
\\Six\\ 



and then taking the limit {5{X, n) — > 0. 



3. Formula (2.13) also holds if 5iX = \\5ix\\t^ and yields 



(i) 



X ~\~ f) $i,X 



in this case. This means the angle a in Fig. 2.15 equals zero, and x, x\ 
is a straight line segment with midpoint y^>. 

The case biX = —\\Six\\t^ is equivalent to (2.14) becoming zero, resp. 
to vanishing of the denominator in (2.15). Geometrically a = tt and 
x^Xi passes through oo. This time multiplying (2.13) with (SiX,v^) 
and using \\t^\\ = 1, one gets = oo in the limit — > 0. 



4. If is a point sphere, its Lie coordinates (A. 2) are also obtained by 



normalizing the eg-component of (2.16) to 1. If cr^ turns out to be a 



plane, Lie coordinates (A.4) are obtained by normalizing (2.16) such 
that the first three components build a unit vector in M 3 . 



3 Cyclidic nets 

The discretization of orthogonal nets as cyclidic nets is closely related to the 
discretizations as circular, conical and principal contact element nets (see 
the introduction). The idea behind cyclidic nets is to construct piecewise 
smooth geometries by gluing cyclidic patches in a differentiable way. For a 
unified treatment we define cyclidic nets in a more abstract way in terms of 
frames at vertices of a circular net, without referring to cyclidic patches. 
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In IR 3 , the realization of a 2D cyclidic net using cyclidic patches is a piece- 
wise smooth C 1 -surface, to be considered as discretization of a curvature 
line parametrized surface (cf. Fig. |1.1| ). A 3D cyclidic net in turn discretizes 
a triply orthogonal coordinate system, by discretizing its coordinate sur- 
faces as 2D cyclidic nets. This corresponds to the classical Dupin theorem, 
which states that the restriction of a triply orthogonal coordinate system 
to any of its coordinate surfaces (by fixing one coordinate) yields a curva- 
ture line parametrization. In the discrete case the corresponding 2D cyclidic 
nets indeed intersect orthogonally along their curvature lines, although the 
definition of a 3D cyclidic net requires orthogonality only at vertices (cf. 
Fig. Oh. 



3.1 Circular nets and definition of cyclidic nets 

We call a function / on Z m a discrete net and write 



where ej is the unit vector of the i-th coordinate direction. Mostly we will 
omit the variable z and denote the functions values by /, fi, fij, etc. 

Definition 3.1 (Circular net). A map x : Z m — > M. N is called a circular net, 
if at each z G Z m and for all pairs 1 < i < j < m the points (x, Xi,Xij, Xj) 
lie on a circle. 

Circular nets are a discretization of smooth orthogonal nets and objects 
of Mobius geometry (see |BS08j for details). For N = 3, the case m = 2 
is a discretization of curvature line parametrized surfaces in IR 3 , where the 
case m = 3 discretizes triply orthogonal coordinate systems. A circular 
net x : 1? — > M. N is determined by its values along three coordinate planes 
intersecting in one point. This is due to the following classical result (see 
for example [BS08J for the proof). 

Theorem 3.2 (Miquel theorem) . Letx,Xi, and Xij (1 < i < j < 3) be seven 
points in~ML N , such that each of the three quadruples (x,X{,Xij,Xj) lies on a 
circle Cij. Define three new circles T^Cij as those passing through the point 
triples (x/;, Xjfc, Xjfc). Generically these new circles intersect in one point: 

X123 = T1C23 n T 2 C13 H T3C12. 

We call an elementary hexahedron of a circular net a spherical cube, since 
all its vertices are contained in a 2-sphere which we call Miquel sphere (cf. 



Cyclidic nets can be defined as circular nets with an additional structure: 
Definition 3.3 (Cyclidic net). A map 



(x, B) : IT 1 R N x { orthonormal m-frames B = (t (1) , . . . , t (m) ) } , m < N 



nf(z) = f(z + ei) = fi(z) 



Fig. 3.1). 




is called an mD cyclidic net, if 
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Figure 3.1: A spherical cube: The eight vertices of a cube are contained 
in a 2-sphere if and only if each face is a circular quadrilateral. Any seven 
vertices determine the eighth one uniquely. 



i) The net x is a circular net for which all quadrilaterals (x 
embedded, 1 < i < j < m. 

ii) The frames B,Bi at neighboring vertices x,Xi are related by 

Bi = (H®oF®)(B) = (F®oH®)(B), 1 < i < m, (3.1) 

where is the reflection ( |2.7[ ) in the bisecting hyperplane (2.6) and 
FW maps (t( 1 ),...,t»,...,t«) to (t^ , . . . , -t (i \ . . . , t^). 



Xijk 




Figure 3.2: Notation for cyclidic nets. 



Extension of circular nets to cyclidic nets. The idea is to choose an 
initial frame Bq at xq, and then evolving this frame according to (3.1 ). This 
is always possible due to 

Proposition 3.4. Given a circular net x : Z m — > M. N with all elementary 
quadrilaterals embedded and an orthonormal m-frame Bq in xq = x(zq), 
there is a unique cyclidic net respecting this data. 
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Proof. As an immediate consequence of Lemma 2.10 the evolution (3.1) is 



consistent, i.e. yields a unique frame at every vertex. □ 

Let k = Y^=i N — i for m < N . As there is a A;-parameter family of 
orthonormal m- frames Bq at xq 6 R , an mD circular net x : Z m — > M. N 
can be extended to a A;-parameter family of cyclidic nets. 



3.2 Geometry of cyclidic nets 

Given the vertices of a cyclidic net, the frames determine a unique cyclidic 



patch for every elementary quadrilateral (cf. Theorem 2.14). The evolution 
of frames guarantees that patches associated to edge- adjacent quads share 
a boundary arc, and supposed all involved patches are non-singular the join 
of patches belonging to the same coordinate directions is C 1 . Nevertheless, 
since the frames are arbitrary, the involved patches may be singular. 

Our main interest is the discretization of smooth geometries, thus in the 
following we assume that all considered patches are non-singular. In par- 
ticular this implies that all vertex quadrilaterals have to be embedded (cf. 



Fig. 2.8), which is already incorporated in the Definition 3.3 of cyclidic nets 



Geometry of 2D cyclidic nets in 



2D cyclidic nets as C 1 -surfaces. Consider a 2D cyclidic net (x,B) in 
M 3 . Identify the orthonormal 2-frames B = (t^ 1 ),^ 2 )) wit h the 3-frames 
(t^\t^ 2 ' ,n), where n = £W x f( 2 ). According to Theorem 2.14, the data 
(x, x\, x±2, X2) and (t^\ , n) at x determines a unique cyclidic patch. 
There are three more circular quads containing the vertex x. Since all quads 
are embedded, for each corresponding patch the vertex frames are related 



by (2.9). Comparison with (3.1) yields that frames fit as shown in Fig. 3.3 



This implies that patches sharing two vertices also share the corresponding 



boundary sphere and boundary arc (see Fig. 2.12 and recall that a circle 
is determined by two points and one tangent). Since the boundary sphere 
is tangent to both patches along the common boundary arc one obtains a 
piecewise smooth surface which is C 1 across the joints. This consideration 
also shows that one should associate circular arcs and boundary spheres to 
the edges of Z 2 . 



Discrete curvature lines. For a 2D cyclidic net the boundary arcs of 
its patches can be treated as discrete curvature lines. These arcs constitute 
a net of C 1 -curves with Z 2 combinatorics on the surface, which intersect 
orthogonally at vertices. The curves are determined by the data (x,B), 
without referring to cyclidic patches. 
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Relation to principal contact element nets. For a 2D cyclidic net 
(x,B) with frames B = (iW,i( 2 )) the normals at patch vertices are given 
by n = x t^ 2 \ Those normals induce a net (x,n) of contact elements. 



As explained in the proof of Theorem 2.14 contact elements (x,n), (xi,rii 



share an oriented sphere (cf. Fig. 2.12). This is the defining property for 



a net of contact elements (x,n) to constitute a so-called principal contact 
element net. Principal contact element nets are a discretization of curvature 
line parametrized surfaces in Lie geometry, and they unify the previous 
discretizations as circular nets in Mobius geometry and as conical nets in 
Laguerre geometry. For details and further references see [BS08]. 

Note that principal contact element nets discretize literally the Lie geo- 
metric characterization of curvature lines, i.e. infinitesimally close contact 
elements along a curvature line share the corresponding principal curvature 
sphere. 

Parallel 2D cyclidic nets. As any C 1 -surface, a 2D cyclidic net (x, B) 
possesses a family of parallel surfaces. In particular, the offset of a 2D 
cyclidic net is again a 2D cyclidic net. The corresponding discrete data is 
obtained by offsetting the circular net x i— > x + en and keeping the frames 
B = (tW.t^.n) (Indeed one obtains another cyclidic net, since the offset 



surface possesses the same reflections (2.7)) 



Ribaucour transformation of 2D cyclidic nets. According to the 
"multidimensional consistency principle" for the discretization of smooth 
geometries, the transformations of cyclidic nets are defined by the same ge- 
ometric relations as the nets itself. For a detailed explanation sec |BS08, 
Chap. 3], where Ribaucour transformations for principal contact element 
nets are introduced in the same spirit. 

Definition 3.5 (Ribaucour transformation of cyclidic nets). Cyclidic nets 
(x,B),(x + ,B + ) : I? ^S 3 x {orthonormal 3-frames B = (t (1) , i (2) , n)} 
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are called Ribaucour transforms of each other, if for all z £ 1? the corre- 
sponding frames B and B + are related by 

B + = F (+)^ B ) = o 

Here H^ + ' denotes the reflection in the perpendicular bisecting plane be- 
tween x and x + , and denotes the orientation change (t^\ t^ 2 \ n) \-+ 
(tW.t^.-n). 

The existence and standard properties including Bianchi's permutability 
property can be derived exactly as for principal contact element nets. Def- 
inition 3.5 means that two 2D cyclidic nets in the relation of Ribaucour 
transformation build a 2-layer 3D cyclidic net. For each pair of correspond- 
ing vertices x and x + there is a sphere tangent to the cyclidic net (x, B) at 



x and tangent to (x + ,B + ) at x + (cf. Fig. 2.12). 



Geometry of 3D cyclidic nets in R 3 

3D cyclidic nets in M 3 are a discretization of triply orthogonal coordinate 
systems. By Dupin's theorem, all 2-dimensional coordinate surfaces of such 
a coordinate system are curvature line parametrized surfaces, i.e. coordi- 
nate surfaces from different families intersect each other orthogonally along 
curvature lines. The discretization as cyclidic nets preserves this property in 
the sense that the 2D layers of a 3D cyclidic net (seen as C 1 -surfaces) inter- 
sect orthogonally along the common previously discussed discrete curvature 



lines (cf. p. 23). Moreover, we will show that the natural discrete families of 
coordinate surfaces of a 3D cyclidic net can be extended to piecewise smooth 
families of 2D cyclidic nets. 

Elementary hexahedra of 3D cyclidic nets. We call an elementary 
hexahedron of a 3D cyclidic net a cyclidic cube. By definition of cyclidic 
nets, the vertices of a cyclidic cube form a spherical cube, i.e. they are 
contained in a 2-sphere. 

Now consider a 3D cyclidic net (x, B) in IR 3 . According to Theorem 2.14 



for each circular quadrilateral (x,Xj,Xjk,Xk) one has a cyclidic patch /W 
with these vertices and the normal t^' in x. Therefore we obtain six cyclidic 
patches as faces of a cyclidic cube: The faces / W, f^\ f^ are all determined 
by the same frame B at x, while the shifted faces are determine d by the 
shifted frames, i.e. is determined by Bi at X{ etc. (cf. Fig. 
Fig. [3751). 



3.4 



and 



Proposition 3.6. Adjacent faces of a cyclidic cube intersect orthogonally 
along their common boundary arc. 



3 CYCLIDIC NETS 



26 




Di- 



B 



B, 



B, 



Figure 3.4: Left: Notation for cyclidic patches associated to elementary 
quadrilaterals of a cyclidic net. Right: Combinatorics of frames of a cyclidic 



net, obtained by the evolution (3.1). 




c 




Figure 3.5: A cyclidic cube. The vertices are contained in a 2-sphere and 
adjacent faces intersect orthogonally. 



Proof. First note that corresponding boundary arcs of adjacent patches co- 
incide: Since all quadrilaterals are embedded, for each face the vertex frames 



are related by (2.9). On the other hand the evolution (3.1) yields frames 



as in Fig. 3.4 For vertex frames of adjacent patches the tangents associ- 
ated to the common edge coincide. Since also the endpoints coincide, the 
corresponding boundary arcs coincide. 

Orthogonal intersection in the common vertices follows immediately from 
the orthogonality of the corresponding frames. Knowing that corresponding 
boundary arcs of adjacent patches coincide, orthogonal intersection along 



those arcs follows from Proposition 2.15 by considering vertex frames of 
subpatches with common endpoints (for example, faces /W and intersect 
orthogonally in each point Xk of the interior of x,Xk, cf. Fig. 3.6). □ 



Proposition 3.7. For a spherical cube there is a 3-parameter family of 
cyclidic cubes with same vertices. These cyclidic cubes are in bijection with 
orthonormal 3-frames at one vertex. 

Proof. This is a special instance of the extension of circular nets to cyclidic 
nets (cf. p. [22b. □ 
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Figure 3.6: Adjacent cyclidic patches associated to different coordinate di- 
rections in a 3D cyclidic net intersect orthogonally along their common 
boundary arc. 



Corollary 3.8. Three cyclidic patches sharing one vertex and boundary arcs 
with combinatorics as in Fig. \3. 7[ can be extended uniquely to a cyclidic cube. 



A 




Figure 3.7: Three faces determine a cyclidic cube. 



Proof. Seven vertices of a spherical cube are already given and by Theo- 
rem 3.2 there is the unique eighth vertex. Since the three patches determine 
the frame at the common vertex, Proposition 3.7 proves the claim. □ 



Discrete triply orthogonal coordinate system. The restriction of a 
3D cyclidic net to a coordinate plane of Z 3 is a 2D cyclidic net, which we see 



as a C 1 -surface (cf. p. 23). We obtain three discrete families of coordinate 
surfaces 



[{x,B)\ Ei{c) | c G Z j , Ei(c) = {( Zl ,z 2 ,z 3 )eZ 3 \zi = c}. 



Three such coordinate surfaces = (x, B) 



vertex x(c\, C2, C3). By Proposition 3.6 any two surfaces m'^fn intersect 



S«(ci)' 



1,2,3, share the 



orthogonally in a common discrete curvature line (cf. p. 23), in particular all 
three surfaces intersect orthogonally in x (cf. Fig. 1.2). These intersection 
curves are considered as coordinate lines of the discrete triply orthogonal 
coordinate system. 
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Extension to piecewise smooth coordinate systems. A non-singular 
cyclidic cube is a cyclidic cube for which opposite faces are disjoint and adja- 
cent faces intersect only along the common boundary arc. If all elementary 
hexahedra of a 3D cyclidic net are non-singular, then the discrete families 
AfW of coordinate surfaces can be extended to piecewise smooth families 
such that the previously described orthogonality property still holds. 

Theorem 3.9. Let fj® be a pair of opposite faces of a non-singular cy- 
clidic cube. For each point x% G x^Xi there is a unique cyclidic patch /W such 
that X{ is a vertex of and /W intersects the four faces fj, f^ k \ fj. 
orthogonally along common curvature lines (cf. Fig. \3.£ty . Moreover, the 
patch /W is non-singular. 




Proof. Given a cyclidic cube and a point 



lines of f^ and /w) one obtains points Xij G xf~x~lj 



G x,Xi, following curvature 
and x ik G Xk~x^ k (cf. 
Fig. 3.8). Take the obtained points Xij,Xik and go along curvature lines of 

the patches f^\ respectively f^\ which gives Xijk,Xikj G a := xJk^Xijk- 
We will prove Xijk = Xikj = Xijki where x^ is the eighth vertex of the 
spherical cube with r em ainin g vertices (x, Xi, Xj,Xk,Xij,Xik, %jk)- According 

this implies the existence of a unique patch 



2.15 



and 



3.7 



to Propositions 
as claimed. 

Denote by Is the Miquel sphere (see Theorem 3.2) determined by the ver- 
tices (x, Xi,Xj,Xk), and by s the Miquel sphere determined by (x, Xi, Xj,Xk). 

Suppose Xijk 7^ Xikj- In this case a C s, as we have three distinct points 
Xjk, x^k, Xikj G a n s , which in turn gives x^k G s. On the other hand x^k 
is also contained in the Miquel sphere s. We have the chain of implications 



Xijk 7~ Xikj O (Z S =t* Xijk G s 



a C s. 
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This proves the claim x^ = x^ = x^ in the generic case a <f_ s. 

The general validity of Xijk = x~ikj can be shown using a continuity argu- 
ment. The crucial observation is that s is fixed by the vertices (x, . . . , x^^) 
of the initial cyclidic cube, whereas the arc a depends on the frame Bjk = 



(cf. Fig. 3.8). The vector t^ is tangent to s if and only 
if fW is tangent to s, due to the evolution (3.1). Now consider a frame B 
at x with £W tangent to s, and as before choose Xj G xT^j- Obviously there 
exists an axis I through x such that rotating B by an arbitrary small angle 
a £ (— e, e) \{0} around I has the effect that the image tw(or) of fM = tW(0) 
does not remain tangent to s. Denote by x^Xi(a) the circular arc between 
x and Xi determined by the vector v^{a) and choose Xi(a) £ x^x~i(a) s.t. 
Xi(0) = Xi. If Xi(a) is continuous in a, the dependent points Xijk(a),xn s: j(a) 
are continuous in a as well. For all angles a £ (— e, e) \ {0} we are in the 
generic case, for which we know Xijk{a) = Xfcjipt). Thus continuity implies 
Xijk{oi) = Xikj(a) in the limit a = 0. 

A cyclidic patch is non-singular if and only if opposite boundary curves 



do not intersect (see Fig£. |2.7 2.8). Since we started with a non-singular 
cyclidic cube, the patch /W i s non-singular. 



□ 



Corollary 3.10. Each pair f^\fj~ of opposite faces of a non-singular 
cyclidic cube can be uniquely extended to a smooth family F^ of cyclidic 
patches such that boundary curves of / W £ F^ are curvature lines of the 
initial faces f^J^J^J^ (cf Fig. 

Moreover one has 



3.8). 



Corollary 3.11. Through any interior point of a non-singular cyclidic cube 
there pass exactly three (pairwise orthogonal) patches /W g F^ l \ i = 1,2,3 



(cf. Corollary 3.10). 



Proof. For any interior point X123 of a non-singular cyclidic cube there exist 
patches /W £ J 7 ^, » = 1, 2,3, such that X123 := f {1) nf^nf^. This follows 
from Corollary |3. 10 using a continuity argument. The patches are unique, 
since patches of a single family F^ are disjoint. Indeed, suppose two patches 
from the same family, say F^ , intersect in a point x* . Then x* is a singular 
point for each patch / 6 J 7 ^ 2 ) U J-^ 3 ** with x* G /, because patches from 
different families intersect along common curvature lines. This contradicts 
Theorem 3.9 since we started with a non-singular cyclidic cube. □ 



Corollary 3.10 implies that the discrete families iWW of coordinate sur- 
faces of a 3D cyclidic net in M 3 can be extended to cont inuo us families Ai^K 
Now taking three 2D cyclidic nets m" G as in Fig. 1 1 . 2\ determines three 
coordinate axes intersecting in xq = mP^ D D m^ 3 \ and according to 



Theorem 3.9 parametrizations of the coordina te a xes induce par ametr iza 
tions of the continuous families Ai^ (cf. Fig. 3.9). By Corollary 
has orthogonal coordinates on an open subset of M 3 . 



3.11 



one 
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Remark on the smoothness. The interior point X123 depends smoothly 
on the coordinates x\, S2, S3. Indeed the points S12, S23, S13 depend smoothly 
on xi,X2,%3, since the faces of a cyclidic cube are smooth surface patches. 
But then also the point S123 as eighth vertex of the spherical cube with 
remaining vertices x, depends smoothly on x\,X2,x,3. 

3.3 Convergence of cyclidic nets to curvature line parame- 
trized surfaces and orthogonal coordinate systems 

Convergence of circular nets to smooth curvature line parametrized surfaces 
and orthogonal coordinate systems was proven in |BMS03j . Starting with 
a smooth curvature line parametrization or orthogonal coordinate system 
/ : D C M m — > R , a one-parameter family f e of circular nets of correspond- 
ing dimension was constructed which converges to / with all derivatves in 
the limit e — > 0. Here the parameter e is the grid size parameter for the 
discretization grid approximating D. Actually more was shown: Not only 
convergence of orthogonal nets as point maps was analysed, but also conver- 
gence of the associated discretization of frames at vertices was proven. In 
the present paper we have given a geometric interpretation of these frames 
as vertex frames of cyclidic patches. Therefore the results of [BMS03J imply 
the following 

Theorem 3.12. i) Given a smooth curvature line parametrized surface 
f : I\ x I2 — > M 3 , where I\,l2 C K are compact intervals, there exists 
a one-parameter family f £ of 2D cyclidic nets converging to f . The 



frames related by (3.1) converge to the corresponding smooth frame of 
f as well. 

ii) Given a smooth triply orthogonal coordinate system f : I\ XI2 x/3 — > IR 3 , 
where I\ , I2, 13 C K are compact intervals, there exists a one-parameter 
family f £ of 3D cyclidic nets converging to f. The frames related by 



(3.1) converge to the corresponding smooth frame of f as well. 
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3.4 Higher dimensional cyclidic nets 



For an analogous geometric interpretation (related to cyclidic patches in 
M 3 ) of mD cyclidic nets in Mr, N > 4, one has to take the following into 
account: Also in M. N the frames at vertices determine circular arcs which 
should serve as boundary arcs of surface patches, but generically the four 
vertices of an elementary quad and the four circular arcs connecting them 
are not contained in an affine 3-space. Nevertheless there exists a unique 
surface patch for the given boundary, which is the image of a cyclidic patch 
in IR 3 under a Mobius transformation: Due to the evolution (3.1 ) the vertices 



and circular arcs associated to one patch are always contained in a 3-sphere 
(cf. Fig. 3.10). The corresponding patch is obtained by identifying this 
sphere with IR via stereographic projection, so the mD case can be reduced 
to the 3D case. 




Figure 3.10: In dimensions greater than 4, the shown curves are generically 
contained in a 3-sphere, which is determined by the vertices and the tangent 
plane at x spanned by t^' and 



3.5 Computer implementation 

The 3D pictures in this work were produced with a Java application based 
on the projective description of Dupin cyclides in Lie geometry. This appli- 
cation is available as a Java Webstart at the webpages of the DFG Research 
Unit "Polyhedral Surfaces'f} The main tool for visualization is the open 
source class library jReality |jG], while the projective model of Lie geom- 
etry is implemented as part of jTEM [ TB| . 

For a given circular net (2D or 3D) in IR 3 one chooses an initial frame 



which is then reflected to all vertices using (3.1 ). As explained in Section 3.2 
one has a unique cyclidic patch for each elementary quadrilateral. In the 
following we will explain how to choose parametrizations of the individual 
cyclidic patches, such that for the induced global parametrization of the 
whole cyclidic net all parameter lines are continuous. 

To get a concrete parametrization of a single patch determined by (X, B) 
one has to chose additional points y-\y^' on boundary circs oc -. x^, ■ x ■ x j • 
Previously we chose y's to be the midpoints (cf. pp. [17] ff.) But even if 



http : / /www .math. tu-berlin. de/geometrie/ps/ software . shtml 
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(jw) touches in the midpoint of general o~v) does not touch 

. So the idea is to choose points yW along two 



s'- ' in the midpoint of i. 



intersecting discrete curvature lines as midpoints of boundary arcs, and then 
to keep track of the intersection points of \ -parameter lines with opposite 



boundary arcs. These new intersection points should then be used in (2.16) 



in order to obtain cr's for adjacent patches. One obtains parametrizations 
of the patches of a 2D cyclidic net, such that all i-parameter lines are 



continuous (cf. Fig. 3.11). 





o yM chosen as midpoint 

n determined via continuity condition 



Figure 3.11: Evolution of the points j/W for initial 
along two intersecting coordinate axes of Z 2 . 



s associated to edges 



To prove continuity of all parameter lines, we verify continuity across 
the common boundary arc of adjacent patches: Consider the induced pa- 
rametrizations over [0, 1] of these common arcs. Since parametrizations are 
quadratic an coincide at three points, they coincide identically. 

In the case of 3D cyclidic nets closeness conditions enter the game which 
have to be satisfied in order to obtain continuous parameter lines. The 
reason is that y's are associated to edges of Z m , and there are different ways 
to determine the same y. But indeed the evolution of y's is consistent, which 
follows from Theorem 3.9 (cf. Fig. 3.8). 



A Lie geometry of oriented spheres 

For the reader's convenience we include a brief overview of the projective 
model of Lie geometry. Probably the most elaborate source is the classical 
book [Bla29j by Blaschke, a modern comprehensive introduction is contained 
in Cecil |Cec92j . A short introduction can be found in jBS08j. 



Subject of Lie Geometry. Lie geometry is the geometry of oriented 
hyperspheres in M. N = M. N U {oo} = S . Points are considered as spheres 
of vanishing radius while oriented hyperplanes are considered as spheres 
of infinite radius. The transformation group of Lie geometry consists of so- 
called Lie transformations. They are bijections on the set of oriented spheres 
and characterized by the preservation of oriented contact (cf. pp. [34f). 
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Mobius geometry as well as Laguerre geometry are subgeometries of Lie 
geometry. Mobius geometry studies properties which are invariant with re- 
spect to Mobius transformations, and accordingly Laguerre geometry stud- 
ies properties which are invariant with respect to Laguerre transformations. 
Mobius transformations are exactly those Lie transformations which pre- 
serve the set of point spheres, and Laguerre transformations are exactly 
those Lie transformations which preserve the set of hyperplanes. 

The Lie quadric. The space of homogeneous coordinates is l^ 1 ' 2 , i.e. 
R N+3 equipped with an inner product (-, •) of signature (N + 1, 2). On the 
standard basis (ei, . . . ,ejy + 3) this product acts as 



Homogeneous coordinates are marked with a hat, except for basis vectors 
which are written in bold face. Points in a projective space we write as 



Definition A.l (Lie quadric). We denote the set of isotropic vectors in 



For a convenient description of the Lie quadric we introduce the vectors 
e o : = 2 ( e N+2 - £n+i) , := - (eAr +2 + ejv+i) , e r := eAr +3 

and replace the basis (ei, . . . , e^v^) by (ei, . . . , ejv, eo, e^, e r ). Note that 
the vectors eo,eoo are isotropic, i.e. eo,eoo G h N+1 ' 2 , and (eo,eoo) = — \. 

Generalized oriented Hyperspheres in M. N are in bijection with points on 
the Lie quadric: 

• Normalized homogeneous coordinates of a proper oriented hypersphere 
s with center c 6 R N and oriented radius r£R \ {0} are 



' l, i=j = l,...,N + l, 
{ei,ej) = < -1, i =j = N + 2,N + 3, 
, 0, 



x = \x\. 





(A.l) 



We follow the convention that positive radii r > are assigned to 
spheres with the inward field of unit normals. 
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• Points x G M w are considered as hyperspheres of vanishing radius, 

x = x + eo + ||a?|| 2 eoo + • e r . (A. 2) 

[v] G Q N+1 > 2 is a point sphere if and only if it is contained the projective 
hyperplane P(e^), i.e. (v,e r ) = 0. 

• The point oo is also a point sphere and has homogeneous coordinates 

oo = eoo. (A. 3) 

[eoo] is the only isotropic point with vanishing eo- and e r -components. 

• Normalized homogeneous coordinates of an oriented hyperplane p with 
normal n G S^ -1 and offset d G R are 

p = n + • eo + 2(feoo + e r . (A. 4) 

[v] G Q N+1 > 2 is a hyperplane if and only if it is contained the projective 
hyperplane P(e^), i.e. (?,e M ) = 0. 

Oriented contact and contact elements. From the Euclidean view- 
point there are the following possible configurations of spheres in oriented 
contact: 

• Two proper hyperspheres are in oriented contact if they are tangent 



with coinciding normals in the point of contact (cf. Fig. A.l) 





Figure A.l: Oriented contact of two proper hyperspheres. 

• A proper hypersphere and a hyperplane are in oriented contact if they 
are tangent with coinciding normals in the touching point. 

• Two hyperplanes are in oriented contact if they are parallel with coin- 
ciding normal, their point of contact is oo. 

• Oriented contact with a point means incidence. A point cannot be in 
oriented contact with any other point. 

Easy to check (see e.g. |BS08j) is the crucial fact 
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Two generalized oriented hyperspheres s± and S2 are in oriented 
contact if and only if (sijSjj) = 0, i.e. if S\,S2 G Q N+1 ' 2 are 
polar with respect to the Lie quadric. 

Given two spheres in oriented contact, there is a unique point of contact. 
At most one of the spheres is a point sphere, thus we have a unique normal 
n at the contact point x. We identify this pair (x, n) E M. N x § N ~ 1 with 
the 1-parameter family of spheres through x with normal n at x, all being 
in oriented contact. Such a sphere pencil is called a contact element (cf. 



Fig. A.2). For x = oo a contact element consists of parallel hyperplanes 



with coinciding orientation. 




^ {x,n) 



Figure A.2: A contact element is composed of all hyperspheres touching in 
a point x with coinciding normal n. 

The projective translation gives 

Definition A.2 (Contact element). A contact element of M N is an isotropic 
line in WP N+1 ' 2 , i.e. a line contained in the Lie quadric Q N+1 ' 2 . We denote 
the set of all these lines by 

£>N+1,2 ._ { IsQtropic Unes in R pJV+l,2| _ 



A quadric of signature (p, q) in a projective space contains only projective 
subspaces U of dimension dimC/ < min(p, q). Thus the Lie quadric contains 
isotropic lines, i.e. contact elements, but no isotropic planes. Equivalently 
there are no three coplanar isotropic lines in Q N+1 > 2 . This means that if 
three oriented spheres are pairwise in oriented contact, they have to belong 
to one contact element. 

Another important implication is the following: For a line L C Q N+1 ' 2 
and a point s € Q N+1,2 ,s L, there is a unique point a £ L such that s 
and a are polar, i.e. there is a unique isotropic line L s B s which intersects 
L (in a). In other words, for each contact element L and a sphere s not 
contained in L, there is a unique sphere a G L such that s and a are in 
oriented contact. 
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Spheres (si, s 2 ) in oriented contact span a contact element inc[si, s 2 ] and 
each contact element contains exactly one point sphere. Homogeneous coor- 
dinates of point spheres are characterized by vanishing e r -component, thus 
the contact point of a contact element L is the unique intersection LnP(e^). 
Solving (a\r 1 + a 2 r 2) e r = and normalizing the eg-component yields 



Lemma A. 3. Normalized homogeneous coordinates (A.2) of the unique 



point sphere contained in a contact element inc[si,S2] are given by 

x = — - — (r 2 ?i - ri? 2 ) . (A. 5) 

r 2 - ri 

In case that one of the spheres is a hyperplane, i.e. the contact element is 
given as inc[s,p], this formula reduces to 

x = s" — rp. 

Lie sphere transformations. Lie sphere transformations of M. N map 
spheres to spheres and are characterized by the fact that they preserve 
contact elements. In the projective model they are described as projec- 
tive transformations of MP^ +1 ' 2 which preserve the Lie quadric and thus 
map isotropic lines to isotropic lines. 

Description of circles. The description of points on a circle in the pro- 
jective model of Lie geometry is essential for us. 

Lemma A. 4. Four points in M. N lie on a circle, if and only if the corre- 
sponding points in Q N+1 ' 2 n P(e^) C MP Ar+1 ' 2 are coplanar. 

The corresponding Mobius geometric claim is a classical result, see for 



example [BS08, Thm. 3.9]. Lemma A. 4 follows immediately since Mobius 
geometry corresponds to the restriction of Lie quadric Q N+1 ' 2 to the hyper- 
plane P(e^) of point spheres. 

We call a continuous family of contact elements sharing one sphere an 
isotropic cone. Such a cone describes a continuous curve on the common 



sphere (or plane) corresponding to its tip. Lemma A. 4 implies that this curve 
is a circle if and only if the isotropic cone is contained in a 3-dimensional 
subspace of MP^ 1 ' 2 . Actually circles on a sphere s are in bijection with 
planes II C P(e^) n pol[s]. These planes in turn are in bijection with 3- 
spaces V C pol[s] since s £ II, i.e. V = inc[s, II]. As each contact element 
(x, n) containing s may be written as inc[s,x], one has {inc[s,x] | x G II} = 
V n Q N+1 ' 2 = C . In particular C is an isotropic cone. This proves 

Proposition A. 5. Circles on a sphere s are in bijection with 3-dimensional 
subspaces V C pol[s] containing s. The contact elements along such a circle 
are the generators of the corresponding isotropic cone C = V n Q Ar+1 ' 2 . 
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